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Thomson scattering of CMBR temperature anisotropies will cause the spectrum of the CMBR 
to differ from blackbody even when one resolves all anisotropies. A formalism for computing 
the anisotropic and inhomogeneous spectral distortions of intensity and polarization is derived in 
terms of Lorentz invariant central moments of the temperature distribution. The formalism is non- 
perturbative, requiring neither small anisotropies nor small metric or matter inhomogeneities; but 
it does assume cold electrons. The low order moments are not coupled to the higher order moments 
allowing one to truncate the equations without any loss of accuracy. This formalism is applied to 
a standard A-CDM cosmology after reionization where the temperature anisotropies are dominated 
by the Doppler effect for the bulk motion of the gas with respect to the CMBR frame. The resul- 
tant spectral distortion is parameterized by u ~ 3 x 10 -8 , where in this case u is observationally 
degenerate with the Sunyaev-Zel'dovich (SZ) y parameter. In comparison the expected thermal SZ 
y-distortion from the hot IGM is expected to be > 30 times larger. However at z > 5 the effect 
described here would have been the dominant source of spectral distortions. The effect could be 
much larger in non-standard cosmologies. 
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I. INTRODUCTION 

The cosmic microwave background radiation (CMBR) 
is observed to have a spectrum extremely close to a black- 
body (a.k.a. thermal or Planckian) spectrum with tem- 
perature T = 2.725 K [H,[l HH. In addition to contam- 
ination from foreground radio and far-IR sources, devia- 
tions from a thermal spectrum is observed in the direc- 
tion of concentrations of hot gas (galaxy clusters) due 
to the thermal Sunyaev-Zel'dovich effect (tSZ) [5J. If one 
had better sensitivity one should see some amount of tSZ 
spectral distortions everywhere since one expect there is 
ionized gas with T g > 10 6 K along every line-of-sight. 

In the Earth frame the brightness temperature of the 
CMBR is observed to vary by several mK in a dipole pat- 
tern; however this would disappear if the observer were 
boosted into the "CMBR frame". After accounting for 
this velocity dipole there is also observed residual pri- 
mary anisotropies at the 10s of /iK level concentrated on 
the sub-degree scales. These anisotropies are expected 
and observed to be very close to blackbody i. e. the tem- 
perature varies with direction on the sky but the spec- 
trum in each direction on the sky is close to a blackbody. 
If one observes the CMBR with a finite beam instrument, 
because of the anisotropies, and because the average of 
a blackbody spectra with different temperature is not a 
blackbody, the observed spectrum will exhibit a spectral 
distortion from blackbody @. This "beam mixing dis- 
tortion" is very similar to, and can be confused with, the 
tSZ distortions mentioned above. However this is not the 
topic of this paper. 



Here the mixing of anisotropic blackbody spectra due 
to scattering of the radiation is examined. These distor- 
tions will occur even with arbitrarily fine angular resolu- 
tion; and thus are called "resolved spectral distortions" 
in contrast to the " unresolved spectral distortions" from 
beam mixing. As we shall see, apart from a small and 
isotropic primordial tSZ effect, the anisotropic blackbody 
approximation is correct to 1st order in the amplitude of 
inhomogeneities in the universe, and the resolved distor- 
tions caused by coupling of scattering and anisotropies 
arise in 2nd order. Just from this consideration one can 
expect the spectral distortion to be very small. 



A. Roadmap 

Here is a brief description of what is included in this 
paper 



£111 At Temperature transform representation. 

£111 Bt Decompose transform into moments. 

£}Al Practical method to compute moments. 

£111 Ct Central moments shown to be Lorentz invariant. 



£111 Dt Moments are coefficients in Fokker-Planck ex- 
pansion. 



'Electronic address: stebbins@fnal.gov 



£111 Et Some properties of expansion. 

£111 Ft Relation to commonly used power law moments. 

£111 Gt tSZ ^/-distortion in terms of moments. 

£jB} Global frames: 3+1 description of space-times. 
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£1111 At Redshifting of photons in terms of global frame 
velocity gradients. 

qiHBt Thomson cross-section as a convolution. 



£1111 Bt Apply convolution to the temperature trans- 
form. 

£1111 Ct Apply convolution to moments. 



£HV At Spectral distortions are driven by temperature 
anisotropics. 



ElIV Bt Cosmological initial conditions and origin of 
spectral distortions. 



ElIV Ct Order in cosmological perturbation theory 
where different terms become non-zero. 



Generalizes all of the above to polarized light. 
qiVDt Relates T to usual AT/T language. 
£HV Et Lowest order equations for small perturbations. 
£HV Ft All moments for small optical depth. 
qiVGt Adds tSZ ^-distortions. 

i fvl Application to expected mean spectral distortion. 

^D} Details for computing mean distortion. 

i jVIt Discussion of results. 

Many of these sections may be skipped as technical de- 
tails, depending on your interests. 

II. REPRESENTATION OF THE SPECTRUM 

Usually one quantifies the flux of photons by a 
brightness [21] I v [c, x, t] where v, x, t and c are respec- 
tively the frequency, spatial position, time, and direction 
in which the photons are traveling (N.B. usually one uses 
n = — c). The quantities we will end up dealing with are 
a highly transformed form of /„. The first step is to 
transform to the dimensionless quantum mechanical oc- 
cupation number n[v, ■ ■ ■] = c I v \ ■ ■]/ (^ttTiv^) where h 
is the Planck constant. Implicit in these definitions is a 
choice of rest-frame. Of course v will be Doppler shifted 
from one rest-frame to another, but the quantity n is 
rest-frame independent (i.e. a Lorentz invariant) once 
one takes into account the Doppler shifting of the argu- 
ment v. 

Other common representations of the spectra is the 
brightness temperature, Tb, defined by Tb = hvf (fee ln[l+ 
l/n] ), which gives the temperature of a blackbody which 
would produce the equivalent brightness at that fre- 
quency Here fee is the Boltzmann constant. If n ^> 1 
this reduce to the Rayleigh- Jeans temperature, Trj = 
hvn/k^. Both Tb and Trj are frame-dependent. 



A. The Temperature Transform 



A blackbody spectrum has frequency dependence 



,bb 



n BB \x] 



1 



(1) 



where T is the temperature. If for a given x, t, c, the 
frequency dependence differs from this functional form 
we say there is a spectral distortion. 

Next define the temperature transform @, 0, H, 0, 
g[ln[T] ■ • •], of n[v, ■ ■ •], defined implicitly by 



dln[T]q[\n[T},' 



BB 



hv 



knT 



(2) 



In §2 it is shown that g[ln[T], • • •] is sufficient to represent 
the spectral distortions that we are interested in here. 



B. Logarithmic and Central Moments 

It is useful to characterize g[ln[T], • • •] by it's logarith- 
mic moments 



d\n[T] ln[T]" g [ln[T],. 



(3) 



It is shown how to compute from n[v, • • •] in ap- 
pendix From these moment we will define other pa- 
rameters beginning with the grayness parameter 



(4) 



If ^(o) then q[ln[T], • • -]/?7(o) integrates to unity and 
can be thought of as a probability distribution function 
(pdf) for ln[T]. Using this pdf one can define the average 



(f[T, ■ ■ ■]) 



^[T]«H /[T , 
»7(0) 



(5) 



so that 77(„) = ?7(o) ( m PY')- Define the mean logarithmic 
temperature as 



T = e < ln I T ]> = e rl( - 1 > /n '° 



(0) 



and the central moments 



l(n) 



In 



E 



(-l) fe n! r/ {k) ( ?7 (1) 



k\ (n - k)\ ?7(o) V»7(0) 



fc=0 



n—k 



(7) 

Note that by definition it( ) = 1 and um) — 0. The most 
important central moment is the 1st non-trivial one, it(2), 
and for this reason we use the special notation 



U = U(2) 



(8) 



and give u the name: width of the temperature distribu- 



tion. A graybody spectrum has u 



(n) 



for n > and 



g > and this is why g got the name "grayness" . 
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C. Lorentz Transformations 

Under a Lorentz transformation, from frame 1 to 2, 
the various quantities are Doppler shifted 



V-2 



T 2 = 



1 + Z\1 
fx 



(9) 



g 2 [ln[T], 



1 + ^12 

ni[(f + z 12 ) v, ■ 
qi [bx[T]+W 



Zl2\ 



The parameters g and M(„) are Lorentz invariant, and this 
is the main reason why it is convenient to parameterize 
the spectral distortion by g and the «(„). 



D. Fokker Planck Expansion 

Another reason why the central moments are useful is 
that the temperature transform of the background radi- 
ation spectra will generically give a very narrow distri- 
bution of temperature, corresponding to a small overall 
spectral distortion. That being the case, to the extent 
that the blackbody spectrum n BB is well represented by 
the Taylor series expansion of it's argument, v/T, then 
the spectral distortion is given by the moments of the 
temperature distribution. A mnemonic representation of 
this is to imagine expanding q as a sum of derivatives of 
Dirac ^-functions centered on some fiducial temperature 
T , i.e. 



q[ln[T],. 



( — 1 ) m 



rn—0 



d m 5 


In 


T 
_ T °_ 




d\n[T 


m 



(10) 



Substituting this into the temperature transform 
one finds, by integrating by parts and assuming 
lim ln [ T ]_> ±00 g = 0, that 



m=0 



ml 



hv 



where 



An^[x] = (-1)' 



t (f n n BB [ 
dln[x] n 



(11) 



(12) 



By computing the moments of eq. (|10[) one can express 
the coefficients d( m \ in terms of g, T, and ui n ) and vice 



'(«) 



^— v mi 
^ n! (m — n)\ 

n—0 v 7 



In 



T 



E 



,1 



ml (n — m)\ 



In 



m=0 



To 
T 



U(n) (13) 

n j 
«(m) 
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(N.B. u ( o) = 1, u w = 0). 




FIG. 1: Shown are the spectral perturbation from the lowest 
order terms in the Fokker-Planck expansion. In particular 
the thick solid curves show defined in eq. (|12p . As 

one goes from black to gray the thick solid curves correspond 
to m — (a graybody distortion), m = 1 (a temperature 
or Doppler distortion), m = 2 (a u distortion), m = 3, and 
m = 4. For reference a {/-distortion An sz [i] = —3An^[x] + 
An' 2 '[i] is shown as a dashed curve. Here x = hv/{k-BTo) 
where To is a temperature which will depend on the context. 
Note that x An is proportional to the change in Rayleigh- 
Jeans temperature, Trj = (hv/k&)n. 



Substituting eq. (fl"4f into eq. (fTTj) and using the Taylor 
series 



BB 



hv 



3=0 J ' 



T 



i J 



BB 



dln[v]- 



hv 



one finds an alternative Fokker Planck series 



n=0 



knT 



kBT _ 
(14) 



(15) 



Eq. (|15p involves the physically defined T rather than 
the arbitrary T in eq. (jlip . the later is probably more 
practically useful since precision spectral measurement 
are often with respect to a reference blackbody. Of course 
if To = T, then dt m -\ = (1 — g)ur m ), the two equations 
are equivalent and dm = 0. 

This formalism is most useful when the temperature 
distribution is narrow enough that the first few terms 
of these series give an adequate approximation to the 
spectrum. A truncation of these series at m < N we 
call an order N Fokker-Planck approximation. Strictly 
speaking a Fokker-Planck approximation corresponds to 
N = 2. 
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E. Fokker-Planck Asymptotes 

It is easy to compute the asymptotic form of the spec- 
tra. A low frequency expansion is 



An^[x} = ~ 



^ 2p (1 - 2p) n x 2 P- 1 



^ (2p)! 



(16) 



where the B p are the Bernoulli numbers with values 



B = 1, Bi 
and Bi n — ■ 
expansion is 



— 5, while for integer n > 1: i?2«+i = 
£fc=o ( 2fc )! ( ( 2 2 "n-fc)+i)! ■ A high-frequency 



An {n) [x] = ^2p n [rx] e 



(17) 



where the p n are order n polynomials defined by p n [x] 



(-l)"e 



x d n 



-e x with limiting values lim^^ooPn^] 



d\n[x\ 

x n and p n [0] — S n ^o- Both expansions have fairly 
good convergence properties. From this one sees that 
lim^o An (m) [x] = 1/x and lim x ^oo An< TO ) [x] = x m e~ x . 
The asymptotes are always positive but An'™' [x] will go 
negative for n > 4. 



F. Power Law Moments 

The moments (T p ) are also useful, for example the 
bolometric brightness is cx (T 4 ) while the number flux of 
photons is cx (T 3 ). Using eq.s (I10I14[) one finds 



(T p ) = T P 



P 



m=0 



in 



00 n 

(m) = (l-g)f»Y?-^ . (18) 



n=0 
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G. Sunyaev Zel'dovich Distortions 

While the tSZ effect is not the focus of this paper, it 
does fit neatly into the formalism just described. The 
Kompane'ets equation which gives the evolution of spec- 
tral distortions due to the tSZ effect, is an N = 2 Fokker- 
Planck approximation to the collisional Boltzmann equa- 
tion with Compton scattering by non-relativistic elec- 
trons. A cloud of hot electrons illuminated by a black- 
body spectrum with temperature To will emit photons 
with a distorted spectrum of the form n « n BB [x] + 



y An 



sz 



Thomson cross-section, n e \ is the space density of free 
electrons, and k B T c is 2/3 of the average kinetic energy 
of the electrons (so T e gives the electron kinetic tem- 
perature for a Maxwellian distribution). The validity of 
eq. HU) is y, t < 1 and k B T c , k B T Q hv < m c c 2 . The lat- 
ter requirements give the condition for the width of the 
scattering kernel in frequency space to be narrow, which 
is the requirement for a Fokker-Planck approximation. 

One expects to find a varying ^/-distortion as one scans 
across the sky, varying according to the density and tem- 
perature of hot gas along the lines-of-sight. From cq. (fTT 
one sees that 



fJ 
f 

U = U( 2 ) 



(20) 



1 + 0[y 2 } 

T (l-3y + O[y 2 }) 

2y + 0[y 2 } 

+ O[y 2 ] n>3. 



From eq. (|18|) one recovers the well known result that 
the change in the number flux of photons is unchanged: 
= (T 3 ) /Tq 3 — 1 = 0, while the increase of bolometric 

brightness is 4r =^ 4 lZf(Lr 1 = %• 

Various authors [HI EH [H, El have examined higher 
order corrections in O e = k B T e / (m e c 2 ), by expanding the 
Boltzmann equation in powers of O c , obtaining expres- 
sions like eq. (fTT) except that the coefficients are pow- 
ers of 8 C . Note that this is only a rearrangement of 
eq.s (|11I15I) and the spectral shapes obtained are linear 
combinations of the ArS m \ For r, hv <^ 1 and a con- 
stant temperature gas, one express the spectral distortion 
by eq. I|15p where the U(„) are given by a power series in 
e beginning with power Q c n l - 

The tSZ effect produces a generic Fokker-Planck dis- 
tortion, however it can be distinguished from other 
sources of spectral distortion, even those described by a 
low order Fokker-Planck expansion. The reason for this 
is that the tSZ distortions are a one parameter (0 e ) set 
of distortions, so that the Fokker-Planck coefficients are 
not independent but are correlated. From eq. (flU)) one 
sees that for a pure y-distortion (©[Gg 1 ]), one expects 
ln[T/To] = — 3u( 2 )/2- Even if it were contaminated by 
primordial anisotropies with U( 2 ) <C ln[T/To], one could 
expect to find the unique tSZ signature: 



Cov 





T ' 




u (2 ),ln 







■-Vax[u (a) ] (21) 



Other sources of distortion will not have this correlation. 



where 



An 



SZr 



" (e*-l) 2 
= -3An (1) fx 



1 



e x - 1 
■An^lx] 



(19) 



which is called a y-distortion. Here the "y parame- 
ter" is given by a line integral through the gas y = 
J drk B (T c — T )/(m e c 2 ) where r = f cn e \UT dt gives the 
Thomson optical depth, ctt = 6.65 x 10~ 25 cm 2 is the 



III. BOLTZMANN EQUATION 

So far we have just discussed mathematical descrip- 
tions of the spectral distortions but not how they arise 
or their physical dynamics. The equation which describes 
the evolution of the distribution of photons in phase space 
(space, time, direction, and frequency) in the presence of 
scattering is the collisional Boltzmann equation. It can 
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be written as a dynamical equation for the occupation 
number n: 



D _dr 



(22) 



where Mi is a convective derivative along a null geodesic, 
C is a dimensionless collision term, and % is the rate of 
increase of scattering optical depth along the geodesic. 
This says that n is conserved along trajectories in phase 
space except for the effects of scattering (which is really 
the definition of scattering). In an arbitrary space-time 
the trajectories are given by the geodesic equation, while 
C and ^ depends on the type of scattering. 
The convective derivative may be written 

D d dx dt dln[v] d 

m = dt + ^t x+ dt & + ^Tdh^]- ( ' 

This equation is purposely written in 3+1, Galilean lan- 
guage, but is applicable in an arbitrary space-time when 
one chooses a global rest-frame (see appendix [B] for the 
specifics). Here 4| give the flow in position space, while 
gives the flow in direction space, i.e. the bending of 
light. These two terms are generally what is known as 
"lensing" , but in this paper the main interest is in the 
4th term, i.e. the flow in frequency space. 



Redshift 



In appendix [B] it is shown that 



d\n[v] dln[l + z] 



dt 



dt 



c ■ (Vv) • c + c • a 



(24) 



where Vv and a are respectively the spatial velocity gra- 
dient and proper acceleration of the global frame we are 
using, while z is the redshift (1 + z is only defined up to 
a multiplicative constant). 

This is not an expression for redshift which may be 
familiar to you, but this non-perturbative expression 
is completely intuitive from well-known, perturbative, 
Newtonian and cosmological results. The 1st term on the 
right-hand-side (rhs) gives the line-of-sight component to 
the velocity gradient of the observer, so this terms can 
be interpreted as a Doppler redshift having to do with 
differences in the observers' velocities. One can interpret 
—a as the gravitational acceleration (g — —a in Newto- 
nian gravity) and the — c • a term gives the gravitational 
redshift as the photons move against or with the "force of 
gravity" . In different frames the amount of Doppler and 
gravitational redshift will differ, but if one considers two 
fixed endpoints the gravitational redshift plus the correc- 
tion for the Lorentz boosts between different frames will 
always agree. 

Eq. (|24p also produces a non-perturbative expression 
for "cosmological" space-times. Consider the metric 



Sa/3 



diag[ 



_ e 2^ e -2V ; g-2V g-2 



</>] 



(l + z c 



(25) 



where z c [rj\ is the cosmological redshift, rj is the confor- 
mal time coordinate, x are the spatial coordinates, while 
4>\x.,rf\ and rj] give arbitrarily large perturbations 
from a flat Friedman-Robertson- Walker (FRW) space- 
time (if c/> = ip — this would be a flat FRW cosmology) . 
In the coordinate frame, the solution of eq. (f24|) is 



1 



Z-2 



1 



1 + z x 1 



\r\ 2 \ 4>i-4>2+ P 2 <M<H-V>) 

-i — re J, 'l 



(26) 



for the redshift between two points 1 and 2 along any 
null geodesic, r[rj\, in this space-time. Here ' = ^ and 
<f>i = <fi[r[r)i], rji\. Note that the logarithm in eq. ([24]) 
means that sum of terms on the rhs translate into a 
product of terms in the expression for 1 + z. Note 
that if (j> = ip = then then the coordinate frame is 
free-falling, so a = 0, and there is only a Doppler red- 
shift. More generally the Doppler term is divided into 

two factors: !""|~ Zc [' )2 j and eJ** 7 ^; the first is known as 
the cosmological redshift given by the velocity gradient 
when <j) = ip = and the 2nd term arises because when 
ip the coordinates expand or contract leading to ad- 
ditional Doppler shifts for observers moving with the co- 
ordinates. The gravitational redshift is also divided into 
two factors, by splitting c • a into a term which a per- 
fect derivative along the null geodesic and yields e^ 1- ^ 2 

and the remainder which yields ei* dl7 ^. Mimicing per- 
turbative cosmology terminology one would call e^ 1- ^ 2 
the Sachs-Wolfe (SW) redshift and e ^ 2 dv (<f>+4>) the in- 
tegrated Sachs-Wolfe (ISW) redshift. These expressions 
are remarkably similar to the well-known result for lin- 
ear perturbations; even though this is completely non- 
perturbative. Of course eq. (|25|) does not include the 
most general perturbations from a flat FRW space-time 
(i.e. no vector or tensor perturbations), and one still has 
to solve for the trajectory r[?7] first. 



B. Thomson Scattering in the Baryon Frame 

Now since we are using a particular global frame in 
which to write the Boltzmann equation, it is most conve- 
nient to use the baryon frame in which the electrons are 
at rest. An important point here is that we are consider- 
ing only the spectral distortion for a cold gas of electrons, 
and specifically ignoring the electron velocity dispersion, 
which will also produce a spectral distortion via the tSZ 
effect. In some, but not all cases, the tSZ effect will dom- 
inate the spectral distortion considered here, but as the 
anisotropy-scattering coupling considered here has usu- 
ally been neglected it is worthwhile to consider it on it's 
own. 

Here we consider not only cold electrons but also the 
limit of low energy photons i.e. hv ^ m c c 2 . In this case 
the cross-section for scattering an unpolarized beam is 
proportional to 



E[c,c'] 



16tt 



(1 



(c ■ c') 2 ) 



(27) 
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Using the the notation 

£*F[a,b r \= f(fc'E[c,c']F[a[c],b[c']]. (28) 



(29) 



the Boltzmann eq. (f2"2"| becomes 
D dr 



—— n = — E*(n — n). 
Dt dt v ' 



Henceforth a "d quantity indicates that it is evaluated at 
direction c'. Note that £-* is a linear convolution operator 
since E*l = 1. 



C. Temperature Transform of Boltzmann Equation 

Substituting eq. ([2]) into eq. (|29|) one finds that the 
Boltzmann equation for g[ln[T], c, x, t] in the baryon 
frame is 



Dt dt w u 



(30) 



where the convective derivative in temperature space is 
D d dc dc dln[l + z] d 



— = 1 V X H V P - 

Dt dt dt x dt 



dt dln[T] 



(31) 



D. Boltzmann Equation of Moments 

When one deal with temperature moments the temper- 
ature dependence is already "removed" so in what follows 
we use the convective derivative 



V 
Vt 



d_ 

dt 



rfx 
~dl 



dc 

~dl 



(32) 



and any effect of redshifting is shifted to the rhs of the 
equation. Since g and are Lorentz invariant there 
will be no redshifting terms for these quantities. With 
this convention if we substitute eq. ([3]) into eq. ([3H)l then 
we find 



Vt 



cflnfl + z] dr ~ . , . . . 

n ~ V(n-l) + -jj- E*W(„) ~ V(n)) ■ (33) 



dt 



This equation tells us that the 7j(„) moments depend only 
on moments with smaller n so that one may without loss 
of accuracy truncate the evolution of the moments at any 
order n. 

What one really wants is the evolution of g, T, the 
central moments ur n y, which one obtains by combining 
eq.s (I4I6I7|) with eq. Q to obtain 



Vg_ 

Vt 

vt 



dt \ J - g 



1-9' 
1 



hi 
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Vt 



u, 



/ » — , n.i 



(n) 



^ ml(n — m)\ 



(n) 

In 



~T'~ 


7i— m \ 


~T 


U '(m)j 



—n In 



U(n-l) - W(n) 



(34) 



Eq.s (|34p . along with the polarized version, eq.s (|C21[) . 
are the main results of this paper. 

As with the r)r n \ this set of coupled equations can be 
truncated at any order N, leading to an accurate rep- 
resentation of all the spectral components An'"' up to 
order TV. Of course a truncation does not give a com- 
plete description of the spectrum. Note that this is not 
a perturbation expansion, rather these are completely 
non-perturbative equations for arbitrary space-times and 
for large spectral distortions. The assumptions are cold 
electrons k&T e -C m e c 2 , soft photons hv -C m e c 2 , only 
Thomson scattering, and unpolarized light. All of these 
assumptions are liable to be real limitations in applica- 
bility. The lack of polarization in these equations was 
done for simplicity, the formulae including polarization 
is given in appendix [Cl In most applications polarization 
will lead only to a small correction and in some cases po- 
larization is completely negligible. These equations are 
tied to the gas frame which may not be most convenient 
for every application. It is relatively simple to translate 
these equations into a different frame: note that g and 
U< n ) are frame invariant while t, c, T and all the dot 
products are not. 



Boltzmann Equation for u 



As we shall see the most important distortions are the 
lowest order ones, i.e. g, T and u — W(2)- The equations 
for g and u are explicit in eq.s (|34[) and here is the explicit 
equation for u = Urn- 



Vu dr ~ 

= — E* 

Vt dt 



(u' -u) + In 



11 

T 



(35) 



One sees that apart from convection, the u distortion is 
sourced by spectral anisotropy, v! — u, and by temper- 
ature anisotropy In [T"/T]. This is illustrative of all the 
moments moments u/ n -\ in that ln[T'/T] can directly pro- 
duce u/ n \ and in that the equation is linear in U( n y For 
n > 2 the scattering term is more complicated and al- 
ways includes non-linear coupling between ln[T"/T] and 
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IV. SPECTRAL DYNAMICS 
A. Temperature and Spectral Anisotropics 

Here the word anisotropy is used to mean any function 
of the n\v, £,•••] at which is zero when n is independent 
of c. This includes functions of q, r)r n -\, g, T, or Ui m \ 
which are zero when there is no c dependence. It is clear 
for eq.s ()30|33|) that the collision term in the Boltzmann 
equation is an anisotropy. One can also see this most 
clearly for eq. ([34]) by grouping the it( m ) dependence of 
[/(„), i.e. decomposing U (n) = £m=o ^( n ,m) where 



u {nAr o 



~T 



(36) 



(n,n)-«(„) — "(n) 

n — m 



u. 



hi 



- U 

~T 
1 



Kn-l) _ u (n-l)) n > 3 



n 



U( m .n) ■ / si 

'ml [n — my. 



In 



l (m) 



n > 3 
€ [2,n-2] 



which is a 



In eq.s (|34|37|) contains terms like In 
temperature anisotropies and terms like g — g', and u'^ — 
Ut n ) which are spectral anisotropies. Note again that if 
the spectral and temperature anisotropies are zero then 
the scattering term is zero. 



B. Initial Conditions and Origin of Spectral 
Distortions 



At very early times electron and atomic scattering is 
sufficient to nearly completely isotropize and thermalize 
the photon distribution to an isotropic blackbody, so the 
initial conditions (ICs) are 



lim t _> {5, In 



1 u (ri>0)} = 



(37) 



So the spectral distortions evolve from zero and the rea- 
son they exist is because of inhomogeneities. Starting 
with these ICs one sees from eq. (|34|) that 



.9 = 



(38) 



Note however that this result does not take into account 
other radiative foreground which can cause g to vary from 
zero. 

The initial growth of temperature anisotropy is given 

by 



din 



lim. 



rrgg. 



T[C! 



din 



l+z[c 2 



1+zlci 



dt dt 
(c 2 - ci) ■ a p + ci • (Vv) ■ ci - c 2 ■ (Vv) ■ c 2 . (39) 



Thus anisotropies in the baryon frame are initially caused 
by acceleration of the gas (a p ^ usually due to pres- 
sure gradients) and/or by anisotropic gas velocity gra- 
dients (shear). Shear will always be a consequence of 
inhomogeneities in the universe and will inevitably lead 
to temperature anisotropies. One sees from eq.s (|34[) that 
anisotropy will lead to time varying, and hence nonzero, 
spectral distortions ur n \. Note however that scattering 
tends to damp existing temperature anisotropies toward 
zero so the amount of anisotropy and associated spectral 
distortions are highly suppressed until scattering turns 
off at recombination. 



C. Perturbative Analysis 

Consider a one parameter family of solutions to the 
full equations-of-motion (EoMs) for matter and gravity 
as well as the initial conditions (ICs). When e = the 
EoMs and ICs give the unperturbed background cosmol- 
ogy, and more generally e gives the amplitude of the 
perturbation from the background solution. One can 
Taylor series the EoMs and ICs about e = and then 
solve them at each order using the lower order solutions. 
The ICs may be stochastic. One can Taylor series in 
e about any quantity which depends on the solution, 
e.g. Q — Q^ eP If the smallest p for which the 

j= is N then one says that Q is a perturbation 
variable of order N and denote this by Q ~ C[iV]. For 
two quantities P and Q, if P — Q ~ 0[N] then one sees 
that P< M > = Q< M > for all M < N. 

The solution of our EoMs, eq. (f3"8"]) . tells us that g ~ 
O[oo] although, as mentioned above, additional radiative 
processes will cause a non-zero g at some order. Since an 
unperturbed cosmology is by definition homogeneous and 
everywhere isotropic one must have 



hi 



{0} 



(V x n){°> = (V 6 n)W = ^ {0> = . (40) 



We know from fllVBI that In 
\n\T'/f] - 0[l\ 



1 + Z[C2,"-] 



It follows "that \n[T'/T] n 
given the ICs of eq. (f3"7) one can see from eq.s (|3~4")l that 



0[1] and hence 
Oln] and 



u 



(«) 
and u 



0[n] for n > 2 (by definition u 



(i) 







(o) = 1 ~ O[0] 
0[oo]). Thus for perturbation the- 
ory at a given order N one need only consider for 
n < N. In linear theory, N = 1, there is no spectral 
distortion, only temperature anisotropy. From the scat- 
tering of anisotropies the lowest order spectral distortion 
is u = U(2) ~ C[2], i.e. spectral distortions only appear 
in 2nd order perturbation theory. 

These conclusions depend on our EoMs which ig- 
nore certain radiative processes. In particular only non- 
relativistic Thomson scattering with hvk^,T c <C m c c 2 
has been assumed. Finite temperature and frequency 
corrections are small, but so are the spectral distor- 
tions. Long before recombination these conditions are 
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violated but the scattering tends to damp spectral dis- 
tortions and I estimate these corrections are relatively 
small although formally they do decrease the order of 
the spectral distortions. Even at O[0] there will be fi- 
nite temperature differences between the electrons and 
photons which will lead to a Oth order isotropic spectral 
distortion from the tSZ effect, i.e. the isotropic part of 
( u (n))c ~ C[0], an effect not included in our EoMs. How- 
ever the amplitude of these terms does decrease rapidly 
with n because it is non-relativistic tSZ. Furthermore 
even (it(2))e is very small. Spectral anisotropies from 
this effect only arise through coupling to temperature 
anisotropies, so the formally anisotropic part, is of higher 
order ui n \ — (ut n \)e ~ 0[1], although this again is a small 
effect and decreases rapidly with n. 

The most important spectral distortion which has been 
ignored, is from tSZ at low z caused by shock-heated gas. 
Shock heating is, arguably, a non-perturbative process, 
but nevertheless may lead to the largest spectral distor- 
tion. 



D. 



AT 

T 



Hubble parameter, the only non-trivial part of Oth order 
Boltzmann equations is 



dhi[T c [t]] 
dt 



-H[t] 



(44) 



Defining the cosmological redshift z c [t] = — 1 + 

e I t dt the solution is the usual temperature- 

redshift relation T c [t] = T c [t ]/(1 + z c [t]). 

The only non-trivial 1st order equation is the usual lin- 
earized Boltzmann equation for AZ. in a frame comoving 
with the baryons. This is what is usually used to calcu- 
late CMBR anisotropies. What is new here only enters 
at ~ 0[2], which is where the lowest order spectral dis- 
tortion arise. These lowest order distortion are solutions 
to the 2nd order equations for u = U(oy. 



dt V dt 
dr<°> 



{0} 



dc 



{0} 



t {2} +(^-Vc) (45) 



E* u'W - + 

dt T T 



One normally denotes a temperature anisotropy by 
AT/T. What this specifically means may vary, but in 
many cases Qf- = (Tb — Tq)/Tq where T is some refer- 
ence temperature and Tb is the brightness temperature 
at a particular frequency. If g = ur n -\ = but T ^ then 
from eq.s ()imi4p the distorted spectrum is 

hv 

I 1 + In — I x = —— . 

e*-l\ T n e*-W k B T 

(41) 

If one Taylor expands the dependence of Tb on ln[T/To] 
about one finds 



+ In 


T ' 


x e x ^ 










e x - 1 / 



AT T b - T 



T 
+ I I 



To 

x 



In 



T 



(42) 





T " 


2 




T " 


3" 






+ o 


In 

















Of course T ~ ^[0] and by the assumed symmetry 
j>{o} _ T c [t] which is independent of x and c. It then 
follows that (AT/T)W = \n[T/T B [t]]W and 



The Oth order gradient operators depends on the co- 
ordinate system one uses in the background cosmol- 
ogy, but if one is perturbing from a flat, CIq = 1, uni- 
verse then one can choose comoving Euclidean coordi- 
nates normalized to z c = so that dc/dt^ = and 
((dx/dt) ■ V x ) {0} = c(l + z c [t])c- V co . Eq. (gSJ) is de- 
rived in the frame comoving with the baryons. In 1st or- 
der perturbation theory, to transform to another frame 
only involves correcting the dipole component of tem- 
perature anisotropy for the Doppler shift caused by the 
velocity of the baryons in the new frame, since angles are 
not changed to 1st order and u is frame independent. 

Thus in an arbitrary frame (i.e. gauge) to 1st order 
make the substitution 



Ar {i} AT {1} r 



{i> 



-(.■•v.: (46) 



where Vb is the baryon velocity in the new frame. The 
lowest order equation for u in the new frame is thus 



In 



{i} 



AT 
~~f~ 



•{1} 



AT {i} 
~~f~ 



(43) 



i.e. for small anisotropies In [T'/T] reduces to the con- 
ventional definition of temperature anisotropy! 



E. Lowest Order Spectral Distortions 



Given homogeneity and isotropy of the background 
solution and using ^ ^ = —H[t] where H is the 



du . A 

— + c(l + z c ) c • V co M 

dr - / . /AT' AT (c - c') -v b 

= — E* \u'-u+[ + jf 

dt \ \ T T c 2 



(47) 



where the order superscripts have been dropped, and r is 
used for the Oth order optical depth, to emphasize that it 
is spatially constant. This equation is sufficient for most 
applications, since the anisotropies are small, the spectral 
distortions are even smaller, and the lowest order spectral 
distortion will be, by far, the largest. 
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F. Single Scattering 

A commonly used approximation where the optical 
depth to scattering is small, for example in the model- 
ing of tSZ effect for clusters of galaxies, is to assume the 
photons undergo at most one scattering, and that the 
light incident on gas is isotropic, unpolarized, and spec- 
trally undistorted. Here allow the incident light to have 
anisotropic T but zero g and U(„) (for n > 1). Under 
these assumption most of the terms in eq.s (|34l) are zero 
and one can express the solution as an integral along the 
photon trajectory: 



V) 





~T'~ 


[ drt*\n 









(48) 



The polarized version of this integral comes by substitut- 
ing g = and < U > ( „ ) = 6 rh0 ^T /2 into eq.s (|C21IC22|) and 



(n) 

integrating to obtain 





~T'~ 


j dr(\jx 






Y 



II I 



(49) 



The trace of this is eq. (|48p . Both of these integrals can 
give a good approximation to all of the moments in some 
situations. 



assumed nearly thermal and isotropic. More generally, 
when the electrons are non-relativistic but isotropic, one 
can use ksT c — ^m c v c 2 , where v 2 is the velocity dis- 
persion of the electrons. For T 7 one wants the average 
temperature of the photons scattered into the beam so 
one should use TV = £*T". Thus one should use 



dy 
dt 



dr fee S*(T C - 
dt m e c 2 



T") 



(51) 



One should expect similarities between the tSZ effect 
and the scattering of anisotropics . This is because un- 
derlying both is non-relativistic Compton scattering i.e. 
Thomson scattering. The relativistic corrections to the 
Thomson cross-section in the center-of-mass frame only 
enters when the center-of-mass kinetic energy KE CU1 ap- 
proaches m e c 2 . For microwave photons this would re- 
quire electron energies of several TeV. For lower energies 
the tSZ effect is just the sum of Thomson scatterings 
with varying Lorentz boosts depending on the motion of 
the electrons. Since g and ui n ) are Lorentz invariant, to 
compute even the relativistic tSZ effect one only needs 
to adjust the formulae for a sum of Lorentz boosts. One 
trivial result is that g — is also fixed point for tSZ ef- 
fect. This result does not depend on the thermality of 
the electron velocity distribution, and holds even when 
the electrons arc relativistic. 



G. Adding Thermal SZ 



AVERAGE u 



Much of the formalism developed in this paper is ap- 
propriate for arbitrarily large spectra distortions but re- 
stricted to cold electrons. The tSZ effect described in 
SQTG] will provide a correction to the collision term of 
the Boltzmann equation due to finite velocity dispersion 
of the electrons. This effect is most well known for the 
case of small spectral distortions and non-relativistic elec- 
tron velocity dispersions 22] in which case one gets a y- 
distortion. In this small distortion limit one can simply 
add the y distortion as a collision term in addition to the 
Thomson scattering term already included. To do this 
one need only note the the relation of y to the central 
moments given in eq. (f20|) . We see that there is no modi- 
fication needed for the Boltzmann equation of g and it( n ) 
for n > 3. The only modifications are for T and u = U(o\, 
which in the baryon frame are then corrected by 



L 1 + 2 J 



Vt 



Vu 
Vt 



dt 
dy 
dt 



(50) 



The validity of this equation is only for small distortions 
so one should restrict oneself to lowest order terms as in 
eq. (|47|) . The usual expression for the total y along a 
trajectory is y — J dt (dr/dt) fee (T c — T 7 )/(m c 2 ) where 
r is the Thomson optical depth used previously, T e is the 
electron temperature, and T T is the photon temperature 



Unlike for anisotropics which have, by definition, zero 
mean, the quantity u must be positive. It is therefore 
interesting to derive an expression for the spatial average 
of u, at a given t i.e. u[t] = («) Xl e- For the lowest order 
dynamics this is straightforward and it is easiest to do so 
in terms of the angular power spectrum of anisotropics 
C/[x, t] at each space-time point defined in the usual way: 

f AT 

a(i, m )[x,t] = / d 2 tYi yl , m) [c]* — [c,x,i] 

1 1 

dlM = 2TTT £ l a (*.-)MI 2 • ( 52 ) 

rn— — I 

Here the y"(; m )[c] are spherical harmonic functions. Tak- 
ing the average of eq. (j47|) . including the tSZ correction 
of eq. ([SO)) , and integrating with initial condition u[0] = 
one finds 



u[t] 
Vsz[t] 



ysz[t] + u v [t] + uatM 
' ,.dr k B (T e - T)e,* 



(53) 



d ^2(|v 7b 



dt 
dr 
dt 3 
dr 
~dt 



m c c 

2\ 



1=2 



2tt 



10 



>l.2 



(Ci, 



10 



Here v 7 b [x, t] is the velocity of the baryons in the CMB 
frame. Here the scattering of the dipole anisotropy in 
has been split from the scattering of the higher I 
harmonics in Hat- In the baryon frame ||v 7 b| 2 /c 2 = 
(21 + l)Ci/(27r) so m7 + «at is the expected mean 
square anisotropy in the baryon frame, not including the 
monopole (I = 0), including the dipole (I = 1), and sub- 
tracting of the quadrapole (1 = 2). One expects that 
the random nature of cosmological inhomogeneities is er- 
godic so one can replace the spatial averages with an 
average over realizations, and use C;'s computed by soft- 
ware like CMBFAST [ij. 

The average u is effected by both scattering of 
anisotropics and the tSZ effect, in contrast since 



£*ln 



T 



o 



one finds 



(T) 6 , x [i]oc(l + z c [t]) e - 3 ^M 



(54) 



(55) 



which is only effected by the tSZ effect and not scattering 
of anisotropies. 



A. Degeneracy of u v and ysz 

These mean spectral distortion for the scattering of 
anisotropies, u7 + mat, cannot be disentangled from the 
tSZ y-distortion because we have no a priori knowledge 
of what the pre-tSZ T was i.e. we cannot use eq. (|55|) to 
determine ysz and then subtract it from u. One way to 
break this degeneracy is to instead consider anisotropies 
in u and T to see what the covariance in these two 
quantities are. For a pure tSZ one expects the relation 
eq. (|21[) . but the scattering effect will decrease the covari- 
ance since scattering does not correlate T with u. The 
primary CMB temperature anisotropies will be a signif- 
icant source of contamination and one might find that 
there is not enough sky to provide good enough statistics 
to measure the covariance accurately enough. 



B. u from Early Times 

At times around and before recombination, z c > 1100, 
the anisotropies are small < 10 -5 but the optical depth 
is very high and the electron and photon temperature is 
nearly in equilibrium so the contribution of these epochs 
to eq. (j53")l is not obvious. However I expect that the 
late-time rather that early-time contribution to will 
dominate. 



C. u from Late Times 

Soon after recombination the optical depth becomes 
very small until the time of reionization, believed to be 



Zrei ~ 10; integrating to a total optical depth of ~ 0.07. 
In standard cosmologies, after reionization by far the 
largest contribution to the anisotropies is the dipole from 
the relative velocity of the baryons and the photons, i.e. 
the i> 7 b term. A large tSZ contribution to is also 
produced at late times as the gas will be adiabatically 
and shock heated as non-linear collapse of structure be- 
gin; and there is also radiative heating as stars, quasars, 
and AGN's turn on. Estimates for the average present- 
day tSZ distortion in a standard A-CDM cosmology are 
ygz[to] « 1 - 2 x lO" 6 [H. 

One can estimate the spectral distortion by late-time 
scattering of temperature anisotropies, by including only 
the dipole anisotropies from velocities, u7 and neglecting 
mat- One should not confuse this effect with the kinetic 
Sunyaev-Zel'dovich (kSZ) effect which also is caused by 
relative velocities of the baryons and photons: the kSZ 
effect produces temperature anisotropies not spectral dis- 
tortions. The late-time velocity contribution to u7 is 
given by 



dr v 2 



dt — 

dt 



7b,rms 



(56) 



where i re i is the time of the beginning of reionization, to 
is today, and M 7 b,rms is the rms baryon velocity wrt to 
the CMBR frame. The calculations is straightforward in 
a standard A-CDM cosmology, and is described in <|D] 
obtaining « 3 x 10~ 8 . There are uncertainties in this 
number but it is clear that in the standard cosmology 
mat <C <C ysz- However this was not always the case: 
the u7 effect is dominated by scattering at z ~ z ro i ~ 
10, while nearly all of the ysz was produced at z < 5 
[15j |. So for z > 5 one expects that ^> ysz- In any 
case all of these numbers are much less than the current 
observational limit u, y < 15 x 10~ 6 0|. 



D. Non-Standard Cosmologies 

If one ventures beyond standard A-CDM cosmologies 
with Gaussian inhomogeneities one can imagine that 
there are regions of the universe where the v 7 b is much 
larger than is observed from our vantage point. Spectral 
distortions, u^, can provide a sensitive probe of regions 
of the universe with large velocities such as might occur 
if there are non-Gaussian voids. This is just what is done 
in ref 01. 



VI. DISCUSSION 

Scattering of temperature anisotropies inevitably lead 
to spectral distortions of the CMBR. These spectral dis- 
tortions have been unknown or ignored to date. They 
are second order in the amplitude of primordial inho- 
mogeneities and very small in standard cosmologies; al- 
though in standard A-CDM cosmologies this mechanism 
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was the dominant source of spectral distortions before 
z ~ 5. 

The main result of this paper is the formalism used to 
compute this effect both in the unpolarized case (main 
text) and for polarized light (appendix). By decompos- 
ing the spectrum into logarithmic central moments one 
obtains a hierarchy of equations where the lower order 
terms do not depend on the higher order ones. One 
can thus truncate the hierarchy without loss of accuracy. 
These are non-perturbative, fully relativistic results for 
arbitrarily large spectral distortions and arbitrarily large 
inhomogeneities. The most limiting assumption is that 
the electrons have small velocity dispersion. Deviations 
from this assumption lead to thermal Sunyaev-Zel'dovich 
effects which are only included in an ad hoc way in this 
paper. 

The spectral distortion in the Oth moment is parame- 
terized by the grayness, g. It is shown that the scatter- 
ing of anisotropies, as well as the tSZ effect, leaves the 
primordial value g = unchanged. The 1st moment is 
parametrized by the mean logarithmic temperature T, 
which provides a global (in frequency space) definition of 
the temperature. 

In standard A-CDM cosmologies it is expected that the 
tSZ effect masks the distortion caused by the scattering of 
anisotropies by more than an order of magnitude, at least 
in the angular averaged spectral distortion. It is possi- 
ble that spectral distortions described here may never be 
measured by looking at anisotropies in the spectral dis- 
tortion. Also this effect may allow one to put limits on 
variants of standard cosmologies. 

The formalism developed here incorporates tempera- 
ture anisotropies, polarization, and spectral distortions 
in a single fairly neat package. This might be useful in 
some pedagogical treatments of the CMBR; especially af- 
ter the tSZ effect is incorporated in a less ad hoc manner. 

A. Future Directions 

Here are some directions for future research to extend 
and apply the results of this paper 

• check whether in A-CDM cosmologies one can re- 
move the tSZ contamination by correlating 
anisotropies in u with anisotropies in T. 

• Look for viable cosmological models where the tem- 
perature anisotropies have larger variation than in 
standard models, leading to larger and perhaps de- 
tectable spectral distortions. One such case has 
been done in ref. [lijj . 

• The formalism used here seems naturally suited to 
apply to the tSZ effect, especially relativistic cor- 
rections. It give a simpler path to expressions for 
these relativistic corrections, especially in the case 
of non-thermal electron distribution functions. 

• Further develop techniques to take the temperature 
transform of real spectra. 



• The central moments for scattering of anisotropies 
as well as the tSZ effect are small, meaning that the 
temperature distribution q[T] is narrowly peaked. 
This will not be true of other sources of contam- 
ination such as synchrotron radiation and free- 
free emission. One can imaging developing a new 
method to filter real spectra, roughly corresponding 
to a notch filter in temperature space, that would 
remove these other contaminants with very high re- 
jection. 



B. Other Features 

Here I list some methods and results of this work, 
which although not the main focus of the paper, some 
readers might find more interesting than the main topic: 

• In SfA] is given a general method for inverting mo- 
ments of a broad class of Laplace-like transforms 
which regularizes singular behavior. The regular- 
ization procedure converts the Laplace-like trans- 
forms into well behaved convolutions which can be 
inverted using Fourier methods. 

• In $C] a transverse tensor representation of the po- 
larization is used , which while equivalent to any 
other representation, leads to extremely simple ex- 
pressions for the Thomson cross-section as well as 
simple evaluation of angular integrals. 

• In $B]a 3+1, global frame, representation of phys- 
ical quantities is developed and used throughout 
the paper. While this is not new and may seem 
less generally covariant, I find it leads to simple 
and very intuitive expressions. 

• In mil Al a non-perturbative expressions for tem- 
perature anisotropies (= redshifts) in cosmological 
space-times w/o scattering is given. These closely 
resemble the linear theory expressions many are al- 
ready familiar with. 

• In eq. ID10I a simple approximation to the linear 
growth of inhomogeneities in a A-CDM cosmology 
which is accurate to better than 1% is given. 
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lim^-ooM 2 ] 



where the ai are non-negative and the (3i are non-positive 
this presents a complication. 

To proceed one can define the differential operator 



A.-II(;£-.)ri(;£-A 



(A3) 



so that D z k\z\ goes to zero rapidly for both large positive 
and large negative argument. On the other hand if k does 
fall off rapidly one can set D z to 1 . 
Define the moments 



Cn 



dz z l D z k[z] 
dx x n D x f[x] 
dyy m f[y] 



(A4) 



where the Ki are known in the sense they can be com- 
puted numerically if not analytically. The and Cm are 
related by 



Cn 



^-^ ml (n — to)! 

m=0 v ' 



E 



-n Cn 



^— ' n\ (m — n)\ 

n=0 v 7 



(A5) 



where the kk can be solved for recursively: 



APPENDIX A: LOGARITHMIC MOMENTS 
FROM BRIGHTNESS 

In this paper the quantities g, T, and u/ n ) are used. 
These are defined in terms of the r]i n -\ in eq.s (|4I6I7[) . 
which in turn are defined in terms of the temperature 
transform (7[ln[T] ,•••]. We will not specify how to go from 
a normal brightness spectrum to q, but we now specify 
how to get the ?7(„) from the spectrum. 

Let us begin with some mathematical preliminaries. 
Suppose we have a convolution of the form 



f[x] 



dyk[x-y] f[y]. 



(Al) 



where the convolution kernel, k, is a known function. The 
relation between the moments of / and the moments of / 
is complicated if k does not go to zero sufficiently rapidly 
for large positive and negative argument. In particular if 



lim. 



«+ 



i=l 



(A2) 




fc-1 

\ - Kk-l 



^ o (k-iy. a 



(A6) 



The inverse exists if and only if kq =/= 0. 

Now let us apply this to the temperature transform, 
eq. @ which is of the form of eq. (|A1[) , when one uses 
x = V = ln[T], f = n and / = q. Thus one finds 

k[z] = (e e — 1) _1 and Cm = V(m)- The function k[z] 
does go to zero for large positive argument so n + = 
but lim 2 _,_ 00 /c[z] = e~ z - | so n_ = 2 with /3j = —1 and 
f3 2 = 0. Thus the differential operator needed to regulate 



divergences is D z 



dz 



and 



D z k[z] 



e z e e * (e z (e e * + 1) - 2 (e e 



1)) 



{e e * - l) 3 



(A7) 



With this form, one can show that for large I the ki 
integral is dominated by large negative argument and 
one can show 



lim/^ooK/ = 



\J dx\n[x] l = ~{-l) l l\ 



(A8) 
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KO = +2 


KO = +2 


(-l) 1 ^! = -0.7820 


/ti = -0.52132 


(-l) 2 |rft2 = +2.0217 


k 2 = -2.42382 


(-l) 3 frc 3 = +0.5424 


k 3 = +6.17299 


(-l) 4 |rK 4 = +1.1669 


k 4 = -7.76751 


(-1) 5 |k 5 = +0.9473 


k 5 = +0.74497 


(-1) 6 |tk 6 = +1.0145 


k 6 = +23.7175 


(-l) 7 ffK 7 = +0.9964 


k 7 = -68.4985 


(-l) 8 !rK 8 = +1.0008 


k 8 = +108.645 


(-1) 9 §k 9 = +0.9999 


k 9 = -55.2136 



TABLE I: Listed are the coefficients 



«(fc) 



and used in 



eq. (| A5|) computed for k < 9 for a temperature transform. 
The Kffc) are computed using the integral in eq. (|A5|l (ana- 
lytically for k — but numerically otherwise) and expressed 
divided by the asymptotic form of eq. jASj. The K(k) are 
computed from the Kim using a backwards substitution algo- 
rithm. 



and this limit is approached fairly rapidly. In table [J the 
first few numerical for the K n and R m are given. 

To go from brightness to r)r n ) for n < N first com- 
pute the occupation number n[v] = c 2 /„/(2/«/), then 
compute the moments 



dv r hv, n d ( 7 dn 

— Ht «- rV 



(A9) 

5|) and from table Q] to 



for n < N. Then use eq. 
compute Cm = Vim)- 



APPENDIX B: GLOBAL FRAMES 

In this paper the occupation number as n[u, c,x, t] is 
used, and the notation indicates a specific spatial coor- 
dinates x, and temporal coordinates t, has been chosen, 
as well as a specific rest-frame in which v and the di- 
rection of travel c is measured throughout space-time. 
This represents a particular choice of a global 4-velocity 
field for the "observer" . Denote the 4-velocity field by 
U a , which is normalized u a g a /su^ — — 1 where g a p is the 
metric. The 4-velocity space-time gradients are canoni- 
cally decomposed into an expansion rate, a shear tensor, 
a vorticity vector, and a proper 4-acceleration which are 
respectively 



u' 
1 

2' 



(Bl) 



Ca/3 = ( + W<5; 7 ~ -0 Ps-y ) P° fi 



On 



Here P a p = g a [j + u a p is the spatial projection tensor 
and e a p-ys is the 4-d Levi-Civita tensor. Thus 



1 



U a ;(3 = ^ Paf3 + CTq/3 + W 7 M e a/ 3j8 - a a Up . (B2) 



Note that P a p is involutive, i.e. P a pPP~ ) 



P a g and a a a are symmetric, and o- a n, 



= P% both 
w 7 , and a a are 
purely spatial, i.e. tr Q ^u /3 = w Q u Q = a Q u Q = 0. 

The tensors P a p and e a p-ys allow one to perform 3-d 
vector analysis in the frame u a , e.g. 

.8 



a-b = a Q P Q/3 6 /3 (axb)-c = a a 6 /3 c 7 



(B3) 



In the geometric optics limit photons will follow null 
geodesies and the frequency observed in a particular 
frame is u[X] oc u a ^j^- where A is an affine parameter. 
Alternately one can define the redshift for that geodesic, 
z[X], such that for all photons following that geodesic 
v oc 1 + z[X]. A frame also induces a time parameteri- 
zation on each geodesic: t[X] — — c j X dXu a ^j^, and one 
will need a convective derivative, i.e. the derivative of 
any quantity along the geodesic, denoted by -Sr. In any 
frame 



D\n[v] _ £>ln[l + z] 



Dt 



Dt 



where 



c ■ a ■ c 



dX 



Cq/3 



dX 



c ■ a ■ c 



<i.r" 
dX 



c • a 



dX ^7<5 rfA 



<2a:T 



dA ™7<5 dA 



da; 15 



(B4) 



(B5) 



Note that if for a frame lj 7 = then the 3-d velocity 
gradient tensor is V • v = |0 1 + a. 



APPENDIX C: POLARIZATION 

The description of the photon distribution in terms of 
n is incomplete because a beam a photons can be po- 
larized. Quantitatively this is not liable to have a large 
effect on the evolution of n but it does have some effect. 
Furthermore to predict the spectral distortions of differ- 
ent polarization modes one will of course need to include 
it in the equations. The equations in this appendix are 
not actually used in the paper but we include them for 
reference. 



1. The Polarization Tensor 

A beam of photons traveling in direction c and space- 
time point" can be described by four Stokes parameters 
as a function of frequency Q w , U u , and V v . Here l v 
is the intensity used above, Q v and JJ V parameterize the 
linear polarization and V v the circular polarization. In 
any particular frame, these four quantities can be used to 

define a 3-d, rank 2 real tensor P called the polarization 

tensor. Here the notation T is used to indicate a 3-d 
tensor, usually of rank 2, which is transverse, i.e. c • 

T [c] = T [c] • c = 0. The simplest such tensor is the 
transverse tensor defined by 



I [c] = I — c 1 



(CI) 
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which is the unique transverse involutive ( I • I = I ) 
rank 2 tensor, and can be used to project into the space 

perpendicular to c, e.g. T = I ■ T • I . 

To construct one needs to define, for each c, two 
direction vectors p'''[c] for i = 1,2 which are transverse 
(c • pW = 0), orthonormal pW • pU) — an d have 
handedness given by p^ 1 -* x p( 2 ) = c. The linear polar- 
ization Stokes parameters Q v and U v are defined up to 
a rotation about the c axis, and they are chosen such 
that electric field oscillations in the ±pW direction cor- 
responds to a Q u > 0, while electric field oscillations in 
the ±(p( 1 ) +p^)/\ /r 2 direction corresponds corresponds 

to U v > 0. In this case P is 



T5>,c,x,<] = £l„ lpW«pW+p( 2 »®p< 2 ) (C2) 



~Q V (pW^pW-p^p^ 



+ \u v (pf^p'npp)®^) 



where ® indicates an outer product. For most purposes 

P provides a sufficient description of the radiation field 
in astronomical applications. The three rotational invari- 
ants are 



Tr[P] 
Tr[c x T 3 *] 



(C3) 



VQ 2 + Ul = \/2Tr[P • P] -Tr[P] 2 + Tr[c x P ] 2 

which are, respectively, the intensity, the circular polar- 
ization, and the amplitude of linear polarization. Unpo- 
larized light has Q v = U v = V v = so 



*P*[c ■■■] = llJc 



I [c] 



(C4) 



which is dimcnsionlcss. The occupation number used 

in the main text is n[v,c,x,t] = Tr[N[i/, c,x, £]]. True 
blackbody radiation is unpolarized and has 



W BB r i 1 



e k B T — 1 



(C6) 



where T is the temperature 



3. Tensor Temperature Transform, Moments, 
Fokker-Planck Expansion 

In analogy with eq.s (I2|i'>l) one can define the tensor 
temperature transform tensor by 



<fln[T] 



QpnpV 



gfc B T _ I 

and the tensor logarithmic moments by 



(C7) 



(n) I' ' ' 



dln[T]ln[T]™ Q [ln[T] 



(C8) 



so that ?7(„) = Tr[H (n )], g = Tr[H (0 )] and T = 

e Tr[H (1) ]/Tr[H Co) ] ( gee eqs gjg])). Finally in analogy 
with eq. JUS]) define 



(C9) 



k=0 



k\(n-k)\ l-g 



so ur n \ = Tr[U („)]; and by definition Tr[ < U > ( )] = 1 and 

and Tx[tf m] = 0. The Fokker-Planck series correspond- 
ing to eq.s (|ll|15H4j) are 



since 



id) 



i(2) 



i(2) 



2. The Occupation Number Tensor 

From P one can construct a quantum mechanical oc- 
cupation number tensor 



N \v, c, x, t] 



2hv 3 



P[i/,c,x,i] (C5) 



"One doesn't mean a literal point, but some finite if small region. 
One needs a finite space-time region dictated by the uncertainty 
principles At Ay < 1, and y<5c|Ax| < c to obtain good angular 
and frequency resolution. The size of astronomical telescopes are 
often dictated by these requirements. 



where 

D ( m ) 
1-5 



• • •] = Y D(m) , r " J AnW 



m=0 



fcBT 



(CIO) 



n=0 



E 



^ n! (m — n)! 



■In 



T 



tf (n) (Cll) 



t? w = E 



— ' m\(n — m)\ 

m=0 y ' 



hi 



To 
T 



D 



(m) 



1-3 



and Tq is an arbitrary reference temperature. 
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4. The Tensor Boltzmann Equation 

In vacuum (i.e. without scattering), in the geometric 

optics limit, P is conserved along geodesies once one 
takes into account redshifting (changes in v) and any ro- 
tation of linear polarization (mixing of Q and U which 
preserves Q 2 + U 2 ). All this can be absorbed into a con- 
vective derivative. Define a "basic" convective derivative 
for transverse tensors, similar to that in eq. 



V_^ _ dx 
Vt dt + dt' y 



|.V,¥(C12) 



dip 

~~dt 



T x c - c x T' 



where dip/dt is the rate of rotation (if any) of the lin- 
ear polarization. This last term does not effect the in- 
tensity or the circular polarization. The full convective 
derivative, D/Dt, includes redshifting, but depends on 
the context. In frequency space 



D V <M[l + z] d 
Dt~vi + dt d\n[v] 



(C13) 



When one includes Thomson scattering off of cold elec- 
trons the collisional Boltzmann equation becomes 



Dt dt 



(C19) 



where r is the Thomson optical depth as defined in mi Gl 

The first term on the rhs, gives the light scattered 

out of the beam, while the 2nd term is the light scattered 
into the beam. Clearly the scattered light is transverse, 
as it must be. It is also clear from this equation that 
the scattering of the antisymmetric part of the polariza- 
tion tensor does not couple to the symmetric part. This 
means that circular polarization evolves independently of 
linear polarization and intensity. 

Note that even if the incoming light is unpolarizcd the 
outgoing light will, in general be linearly polarized, i.e. 
scattered unpolarized light will be polarized. The ap- 
proximation used in the main text is to use unpolarized 

light in scattering term N [u, £,-••] = \n\y, £,-••] *T*[c] 
and then take the trace to get scattered occupation num- 
ber, so that 



in temperature space 



Dn 
~Dt 



dr 3 

~dt lOTT 



d 2 c' (1 + (c • c') 2 ) (n' - n) (C20) 



D V d\n[l + z] d 
Dt,~Vt + dt <91n[T 

for the logarithmic moment tensor 

dln[l + z 



^ H («)-^ H («)- 

for ln[f ] 



dt 



^ln[f] = ^ln[T]-«±il = ^ln 
Dt L 1 Vt 1 1 dt Vt 



(n-l) 



T 



l + z 



(C14) 



(C15) 



(C16) 



while D/Dt = V/Vt for g and t? (ll) . 

Non-convective effects include scattering, refraction, 
dispersion, etc.; but in cosmological application the most 
important effect is (Thomson) scattering off of free elec- 
trons. Before proceeding with Thomson scattering define 
the angular average operator 

(F[a, b'})-, = i- J d 2 c' F[a[c],b[c']] ■ (C17) 

Note that the convolution operator of eq. (|28|) is 

±*F[a,tf] = ^F[a,b'} 1 + {t 2 t,)2 y r (C18) 

= UF[a,bf]Tr{T .T' - *T] 



which is where comes from. 



which is the origin of eq. (|27j) . This is correct if there is 
only one scattering of initially unpolarized light, but not 
for multiple scattering. In many applications the single 
scattering approximation is good. In any case one knows 
that the CMBR is at most about 10% linearly polarized 
so the approximate equation is not liable to lead to large 
errors. 

From eq. (|C19[) one can derive the Boltzmann equation 
for g, ln[T] and U ( n) : 



Dln[l - g] 



Vt 



dr_ 
~dt 



1 



-Tr 



I A 



(0,0) 



I 



V 
Vt' 



-In 



T 


= ^Tr 






.1 + *. 


2 dt 



T.(A ( i i0) +A 



(o,i) J 



^ W - -~ Tr 
Vt {n > 2 dt \ 



T • A (0 ^ 0) • I 



(n) 



-n Tr 



(0,1) 



x ^ n! 
^-^ m\(n — to)! 

m=0 v ' 



I A 



(n— r 




where 



If one sets U / 





T'~ 






~T 



(0 



(C22) 



(„) = ui n ) I /2 on the rhs and take the 
trace one will recover eq. (1341) . 
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Neither Gray nor Circular 



where 



One sees from eq.s (|C21j) that if TJ („) is initially sym- 
metric it will remain so. So under the assumption of 
Thomson scattering by cold electrons one sees that no 
circular polarization will develop. Furthermore one sees 

that if one sets g = 1 and U ( ) = 1/2, which is the 
expected thermal initial conditions, then 

^ = I 

vt vt 



-tJ 



(0) 







(C23) 



So as with the unpolarized equations there are no gray- 
ness terms; and the spectrum of the polarization terms 
are given by derivatives of a blackbody, not a blackbody 
itself. Thus one can simply set the circular polarization 
and grayness to zero. Of course this may not remain true 
if other radiative processes are included. 



b. Perturbative Analysis 



Note that the perturbative analysis of § IIV CI is un- 
changed from the unpolarized version of the equations, 
namely 



O[oo] 
0[1] 



~ u (0) 




~T'~ 


111 






~T 





1[n] n > 1 



(C24) 



B 



(k,i) 



= ((l-«/)m[T'] fe tj 



(') 



(C26) 



The important point is that ~B(k,i) is independent of c 
and depends only x and t. One never has to store or 
compute functions of both c and c'. For each space-time 
point, to compute all the B( ifc )'s scales like One 
can then compute the A-(i t k) locally at each c element, 
using eq. (|C25|) . which scales like Then computing 
the scattering terms in eq. (|C2ip is also local in c space, 
scaling as 4fin- Finally the convective derivative is quasi- 
local in x, t, c-space, involving only adjacent points. For 
large # m the largest scaling is to compute eq. (IC25I) . 



5. Gravitational Field of the Spectral Distortions 

An additional complication in truncating the moment 
distribution is that the stress-energy of the photons does 

not involve only the lowest order moments, TJ („). So 
even though in the radiative transfer sector the moment 
truncation is exact and non-perturbative, a truncation of 
the spectral distortion can lead to only approximate in- 
formation about the gravitational field of the photons. In 
many cases however the gravitational field of the photons, 
let alone the spectral distortions is completely negligible. 



APPENDIX D: FORMULAE FOR LATE-TIME 
SPECTRAL DISTORTION 



c. Numerical Scaling 



One could imagine numerically simulating eq.s (IC21|) 
over some space-time volume. Without scattering the 
full radiative transfer computation scales like the num- 
ber of space-time points (#x#t) times the number of 
frequencies times the number of angular resolution 

elements In general scattering can increase this by 

additional factors of However by using the mo- 

ment decomposition one can reduce to a relatively 
small number of moments # m , e.g. # m = 1 if one only 
wants to simulate T. Within the context of Thomson 
scattering by cold electrons one does this without ap- 
proximation. Furthermore the scattering term only adds 
an additional factor of #^ . Thus the total computational 
scaling is #^#t#e#m- O 1 course full relativistic radia- 
tive transfer requires the timestep to be less than the 
light crossing time of the spatial resolution element. For 
non-relativistic flows various quasi-static approximations 
can be used to reduce the number of timesteps. 

To achieve the scaling one can use the relation 



Hfc,i)[c, • • -J = 



1 



(-l)'fc! 
1-9^ »!(*-<)! 



E 



ln[f [c, • • -]fB 



(fe-i,0 I 



(C25) 



In practice one will find it convenient to express the 
eq. (|56p as an integral over cosmological redshift, z, 
rather than t, and one can use the relationship 



dln[\ + z] 
~dt 



-H[z] 



(Dl) 



where H[z] is the Hubble parameter. Here the common 
H = H[0] and h = H /{100 km/s/Mpc) is used. 

To compute the late-time peculiar velocity dispersion 
of the baryons approximate by using the linear theory 
peculiar velocity dispersion given by 



^rms,lin 1^1 



00 



H[z] f[z] D[z] 
H 1 + z 



(D2) 



where z is now the cosmological redshift, D[z] is the lin- 
ear theory growth factor, f[z] = — (1 + z)D'[z]/D[z] w 
^ml 2 ] ' 6 ! and Aj 2 in [fc, z] is the linear theory dimension- 



less power spectrum (see [l7| §16.2) 
approximately the solution to 



D' 



q[z\ 



2(i + z y 



For lo-z, D[z] is 



D[z) = (D3) 
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with boundary conditions D[0] = 1 and D[oo] = 0. Here 
^m[z] = STrGp m [z]/3H[z] 2 is the matter density parame- 
ter and q[z] = — 1 + (1 + z) H'[z]/H[z] is the deceleration 
parameter. In terms of a sum over the different compo- 
nents, c, 



H[z] = H J2 ^co(l + z) 3 ( 1+ ^ 



(D4) 



where fi c o gives the present density parameter of each 
component and w c the equation-of-state (p/(pc 2 )). In 
the concordance (flat A-CDM) model c sums over A (for 
a cosmological constant), m for matter (baryons + cold 
dark matter), r for radiation (photons and light neutri- 
nos). In the assumed flat cosmology fi m o+f2 r o+^Ao = 1. 
The equations of state are w\ = —1, w m = 0, and 
w T = i. 

The power spectrum is approximated by 



Af in [fc,0]=Afc" +d T B BKS 



T = Q ro0 exp 
Tbbks [q] = ■ 



1>0 




(D5) 



1 o 

\ i 'mO 

ln[l + 2.34 q] 



9 Idn 4 / l+3-89 9 +(16.1q)2 
z.oigy +(5.46<j)3 + (6.71 9 ) 4 



WWN =3 
4 = J y Agjfc, 0] W han [k 8 h-'Mpc} 2 . 

where Tbbks is taken from [l8| . The last equation serves 
to normalize A in terms of the observational parameter 
08 < 

Finally one needs the optical depth. The z = electron 
density is n o = 3iJo^bo/(87rGmH) (1 — I^Kc) where mn 
is the mass of the hydrogen atom and Yhc is the helium 
mass fraction. Thus 



df 

— = T H (l + zf Xc [z} 



(D6) 



3 cHo<jt 
87rGmH 



1 - - 2 Y He 



where Xc is the fraction of electrons which are ionized. 
Here is assumed instantaneous reionization at redshift 
z re ; which is implicitly defined by 



where uq = |ro 0q. 

To get numbers out the following numerical values for 
the cosmological parameter: h = 0.73, O m o = 0.24, 
n h0 = 0.0416, y Hc = 0.24, CTg: = 0.756, and r obs = 0.9 are 
chosen based on ref.s ([l|l[2(|). Thus r = 1.85 x 10~ 3 , 
z rci = 10.3, fa = 2.7 x 10" 3 , m = 9.1 x 10" 9 , and finally 



= 3.4 x 10" 



(D9) 



This is very small! 

To see how uncertainties in cosmological parameters 
can change this number, let us first note that the growth 
function is well approximated by 



D[z] 



£1. 



(D10) 



for z < 20. Next note that the integrals of eq.s (|D7|D8j) 
are dominated by z ~ z rc j where the cosmological con- 
stant is negligible and ^l m [z] w f[z] f» 1 while H[z] w 
^(l + zf/^oso 



Tobs « 3 T Jl m0 /2 (1 + Z rci ) 3/2 



(Dll) 



u v » -r /3 2 ^ 8 (Vm~-l) 



or eliminating z ro j one finds 



1 



0.08 



;Po T "m0 




(D12) 



One sees a relatively weak dependence on £2 m o and r b s ; 
while r cx H^qHq is fairly well constrained. By far the 
largest uncertainty is through /3q oc erf. There is con- 
siderably debate as to the value of erg. There are also 
non-linear corrections to the rms velocity which have not 
been included. 

Ref. [HI have given a heuristic formula for the late- 
time tSZ effect which includes only heating from gravi- 
tational collapse is 



r obs = r / r " dz (1 + zf ^- (D7) 
Jo n l z \ 

in terms of the observed optical depth, r b s - Thus one 
finds that 

u v =u f" 1 dz^{f[z]D[z]f (D8) 
Jo H o 



ysz[to] = 1.5 x 10 



^8 



4.1-fi„ 



mO 



1.26-o-g 



,0.756/ VO-24, 

(D13) 

For the choice of cosmological parameters given above 
one finds that the tSZ y-distortion is ~ 75 times larger 
than that from scattering of anisotropies. 



